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PREFACE 
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1.0  INTRODUCTION 


One  of  the  most  widely  used  methods  for  calculating  boundary -layer 
flows  is  that  developed  originally  by  Patankar  and  Spalding  (Ref.  1). 

This  has  been  largely  due  to  the  versatility  and  ease  of  application  of  the 
method.  At  the  Arnold  Engineering  Development  Center,  the  method  of 
Patankar  and  Spalding  has  been  applied  to  a  wide  variety  of  problems, 
ranging  from  laminar  incompressible  flows  to  turbulent  supersonic  flows 
(e.g..  Refs.  2  and  3).  As  originally  conceived,  this  method  was 
capable  of  treating  both  laminar  and  turbulent  boundary  layers  and  in¬ 
cluded  the  transverse  curvature  effect  for  both  internal  (nozzle  or  channel) 
and  external  axisymmetric  flows.  Each  user  of  the  method  of  Patankar 
and  Spalding  has,  of  course,  styled  it  to  his  particular  taste,  and  a  num¬ 
ber  of  significant  modifications  have  been  made.  These  include  the  addi¬ 
tion  of  equilibrium  dissociating  air  thermodynamic  and  transport  proper¬ 
ties,  the  treatment  of  entropy  layer  swallowing  in  computing  supersonic 
and  hypersonic  flow  over  blunt-nosed  bodies,  and  extension  of  the  basic 
method  to  treat  the  three-dimensional  boundary -layer  flow  in  the  wind¬ 
ward  symmetry  plane  of  an  axisymmetric  body  at  angle  of  attack 
(Refs.  4  and  5). 

The  purpose  of  this  report  is  to  present  a  recently  developed  method 
for  computing  boundary-layer  flows.  Although  the  method  to  be  presented 
is  derived  from  the  work  of  Patankar  and  Spalding,  there  are  a  number 
of  significant  deviations  from  and  extensions  to  their  work.  For  this 
reason,  much  of  the  success  of  the  present  method  must  be  attributed  to 
the  previous  efforts  of  Patankar  and  Spalding,  while  any  blame  for  short¬ 
comings  of  the  modified  computation  must  be  borne  by  the  present  author. 

One  of  the  primary  contributions  of  Patankar  and  Spalding  was  the 
introduction  of  the  u  coordinate  for  the  direction  normal  to  the  body  sur¬ 
face.  This  is  a  normalized  von  Mises  coordinate  which  varies  from  zero 
at  the  body  surface  to  unity  at  the  outer  edge  of  the  boundary  layer.  The 
simple  u  coordinate  has,  however,  certain  shortcomings:  one  is  the 
singularity  of  certain  derivatives  with  respect  to  u  when  evaluated  at  the 
body  surface,  and  a  second  is  the  necessity  for  predicting  the  rate  of 
change  of  the  stream  function  along  the  outer  edge  of  the  boundary  layer 
as  the  flow  proceeds.  The  first  of  these  can  be  eliminated,  as  has  been 
done  by  Landis  (Ref.  6),  by  adopting  a  so-called  td2  transformation,  and 
this  approach  has  been  used  in  the  present  work.  Not  only  does  the  use 
of  the  cd2  transformation  eliminate  the  aesthetically  distasteful  singularity 
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mentioned  above,  it  also  yields  more  accurate  numerical  results  accord¬ 
ing  to  Denny  and  Landis  (Ref.  7).  Unfortunately,  the  problem  of  deter¬ 
mining  the  growth  rate  of  the  stream  function  has  no  patent  solution,  but 
through  careful  monitoring  and  controls  it  can  be  adequately  determined. 

A  second  difference  between  the  present  work  and  that  of  Ref.  1  is 
that  an  option  for  considering  the  effects  of  longitudinal  curvature  and 
the  variation  of  the  static  pressure  in  the  direction  normal  to  the  body 
surface  has  been  included  in  the  method  presented  herein.  The  use  of  a 
difference  scheme  derived  from  Taylor  series  truncation  is  another 
major  difference  between  the  method  described  in  this  report  and  that 
presented  in  Ref.  1.  Because  of  the  coupling  of  the  normal  and  longi¬ 
tudinal  pressure  gradients,  and  because  of  the  numerical  scheme  used 
for  turbulent  flows,  the  resulting  finite -difference  solution  method  is  an 
iterative  one,  in  contrast  to  the  non-iterative  method  of  Ref.  1. 

The  next  section  presents  the  governing  equations  and  the  finite - 
difference  method  which  is  employed  in  solving  them,  and  in  the  succeed¬ 
ing  section  are  given  the  results  of  a  variety  of  test  calculations  which 
have  been  performed,  including  internal  and  external  flows  and  flows  with 
significant  normal  pressure  gradients. 


2.0  GOVERNING  EQUATIONS  AND  METHOD  OF  SOLUTION 


In  this  section  are  presented  the  governing  equations  for  compressi¬ 
ble  boundary -layer  flow  and  the  finite-difference  method  used  to  solve 
these  equations.  The  boundary -layer  equations  include  transverse  curva 
ture,  longitudinal  curvature,  and  normal  pressure  gradient  effects.  The 
equations  are  given  first  in  x-y  surface  coordinates,  and  later  the  trans¬ 
formation  to  x-w  coordinates  (using  the  u2  transformation)  is  given. 

In  this  investigation,  the  stagnation  point  form  of  the  boundary -layer 
equations  were  derived  and  solved.  To  the  author's  knowledge  this  has 
not  previously  been  done  using  the  to  or  to2  transformation,  and  for  that 
reason,  the  stagnation  point  form  of  the  boundary -layer  equations  using 
the  to2  transformation  is  also  presented  in  this  report. 

The  results  which  are  given  in  the  next  section  are  all  for  a  therm¬ 
ally  and  calorically  perfect  gas;  however,  the  method  which  has  been 
developed  is  also  capable  of  treating  the  flow  of  equilibrium  dissociating 
air.  This  has  been  achieved  by  using  the  thermodynamic  and  transport 
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property  correlations  of  Cohen  (Ref.  8)  which  have  been  faired  into 
perfect-gas,  constant -Prandtl  number,  Sutherland -viscosity  relations 
at  lower  temperatures.  For  the  case  of  turbulent  flows,  a  two-layer 
eddy -transport  property  formulation  was  used;  a  streamwise  inter- 
mittency  factor  was  used  to  treat  transition  from  laminar  to  turbulent 
flow. 


2.1  BOUNDARY-LAYER  EQUATIONS  AND  TRANSPORT  LAWS 

The  boundary -layer  equations  which  have  been  used  are  essentially 
those  given  by  Maslen  (Ref.  9),  except  that  the  Prandtl  number  has  not 
been  assumed  to  be  constant  and  is,  therefore,  included  within  the  y  de¬ 
rivative  in  the  energy  equation.  When  both  longitudinal  and  transverse 
curvature  terms  are  retained,  the  boundary-layer  equations  have  the 
following  form: 

continuity : 


(pufj)  +  [(1  +  Ky)pvrj]  =  0 


(1) 


x- momentum: 


<9  u 

P*T*  + 


pv(l 


.  du 
+  K  v)  -r- 
dy 


+  pUVK  =  - 


(2) 


y -momentum: 


dp 

dy 


=  Kp  11 


2 


(3) 


energy: 


«9H 

Pu^  + 


P^I+KV)^-  = 


*d(pu2)  d_ 

dy  +  rj  dy 


n  \ 


Pr-  1 
Pr 


prKl  +  Ky) 


d(n2/2) 

dy  _ 


(4) 


These  equations  are  similar  to  those  used  by  Back  (Ref.  10),  although 
he  treated  some  of  the  curvature  terms  differently.  In  these  equations 
u  and  v  are  velocity  components  parallel  and  normal  to  the  surface  (in 
the  x  and  y  directions),  and  k  is  the  surface  curvature  (see  Fig.  1). 

The  k  is  the  reciprocal  of  the  longitudinal  radius  of  curvature  R,  except 
that  it  is  counted  negative  for  concave  surfaces  and  positive  for  convex 
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surfaces.  The  power  j  is  +1  for  axisymmetric  flow  and  0  for  two- 
dimensional  flow.  The  r  and  y  are  related  by 

r  =  rw  +  *y  cos  <*>  ^ 

where  rw  is  the  transverse  radius  of  the  body,  0  is  the  body  slope,  and 
i  is  an  index  which  is  +1  for  external  flow,  -1  for  internal  flow,  and  0  if 
transverse  curvature  is  to  be  neglected. 


Figure  1.  Typical  configuration  showing  coordinates  and  velocity  components. 


Boundary  conditions  (at  the  surface,  subscript  w,  and  outer  edge, 
subscript  e)  for  Eqs.  (1)  through  (4)  are  as  follows: 

y  =  0:  v  =  vw 

u  =  0 

H  =  II  w  or  dll  /dy  =  0 

v  =  ye:  u  =  ue 

P  =  Pc 
II  =  H 

oc 

Typically,  the  distribution  of  static  pressure  along  the  body  surface,  as 
computed  by  methods  for  inviscid  flow,  is  known.  Assuming  then  that 
Eq.  (3)  is  valid  in  the  inviscid  flow  region  near  the  body  surface,  pe  was 
determined  by  numerically  integrating  Eq.  (3)  from  the  body  surface  to 
the  value  of  the  stream  function  at  the  outer  edge  of  the  boundary  layer. 
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treating  the  flow  as  inviscid.  This  was  done  while  simultaneously  de¬ 
termining  the  variation  of  u,  p,  etc.  in  this  region  by  conventional 
methods,  including  the  treatment  of  entropy  swallowing  on  blunt  bodies 
in  supersonic  flow  by  the  method  described  in  Ref.  4.  This  also  re¬ 
sulted  in  the  determination  of  ue. 

For  laminar  (molecular)  thermodynamic  and  transport  properties, 
constant  specific  heats  and  a  Sutherland  viscosity  law  were  used  for 
perfect -gas  calculations  and  Cohen's  correlations  (Ref.  8)  mentioned 
previously  were  used  to  obtain  properties  of  equilibrium  dissociating 
air. 


In  order  to  treat  transitional  and  turbulent  boundary -layer  flow,  a 
two-layer  eddy -transport  property  hypothesis  and  a  streamwise  inter- 
mittency  factor,  7,  were  employed.  The  independent  variables  in 
Eqs.  (1)  through  (4)  were  treated  as  time -mean  quantities,  and  the 
transport  properties  were  treated  as  follows: 


where 


m  M  +  ypi 


=  pk 


2  2 


(7) 

(8) 


(9) 


^Pouter 


>2  2 
PA  >m 


(10) 


Equation  (9)  was  used  out  to  the  point  where  (dinner  =  ^t^outer*  anc* 
Eq.  (10)  was  used  beyond  that.  Equations  (9)  and  (10)  are  those  recom¬ 
mended  in  Ref.  1,  and  the  values  used  for  A,  A*,  and  k  were  0.09, 
26.0,  and  0.435,  respectively,  which  are  essentially  those  given  in 
Ref.  1.  The  ym  is  the  value  of  y  where  u/ue  =  0.  99. 


The  intermittency  factor  determined  by  Dhawan  and  Narasimha 
(Ref.  11)  was  used;  this  has  the  form 


(ID 
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where  xa  is  the  value  of  x  at  the  beginning  of  transition  and  xb  the  value 
at  the  end  of  transition. 


2.2  TRANSFORMATION  OF  THE  GOVERNING  EQUATIONS 

In  order  to  transform  the  governing  equations,  a  stream  function 
was  defined  which  satisfies  the  continuity  equation,  and  then  a  variable 
lj  was  introduced,  where  u2  is  a  normalized  form  of  the  stream  function. 
The  stream  function  ip  is  defined  by 

(fl)  =  -  U  +  K>'W 


The  variable  u  is  defined  by 


Cl) 


2 


(12) 


(13) 


These  were  used  to  obtain  the  transformation  operators 


W 

_  /dt  \  ( 1  +  Ky)p  v  fj  /d(  \  dx 

dx/y  \dco/xZ ~  \chj)  x 


_  p  u  rj  /<9f\ 

\dy/x  2<u(^e  -  i//w)  \do})x 


2  w(0e-^J 


(14) 

(15) 


Equations  (14)  and  (15)  can  be  combined  to  treat  the  convective  terms 
of  Eqs.  (2)  and  (4),  yielding 


pu 


y 


+  (1  +«y)pv^|i^ 


(16) 


Application  of  Eqs.  (14),  (15),  and  (16)  to  Eqs.  (2)  and  (4)  yields  the 
following  general  form 
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where,  for  f  =  u 


b 


u 


C 


u 


d 


u 


puu 


( 1  +  *y) 


2tu  PEI" 


»2(«L»r  - 


p  u  PEI 


(dp\  (l+Ky)vrJ  KurJ  ( dp\ 

+  \fr)x  2(0  0  PEI  ■  2(0  PEI  \dco/x 


with 


and 


d^e 

“d7 


=  -rJ 


r  e  m 


e 


d«Aw 

dx 


rJw  mw 


(17) 


(18) 


(19) 


PEI  = 


w 


w 


For  the  energy  equation,  where  f  =  H, 

a.,  =  a 
rl  u 


li 
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In  addition  to  solving  equations  of  the  form  of  Eq.  (17)  for  u  and  H, 
Eq.  (3)  was  integrated  from  y  =  ye  to  y  =  0  to  obtain  the  variation  of  p 
through  the  boundary  layer,  and  the  definition  of  the  stream  function 
was  used  to  obtain  the  v  component  of  velocity.  To  obtain  y  from  u,  the 
following  equation  was  integrated  at  a  given  value  of  x  from  w  =  0  to 

to  =  1 


2t>  PEI  , 

y  =  /  - r 

0  Pur)  (21) 

The  quantities  m''  and  m^'  introduced  in  Eqs.  (19)  were  required 
in  order  to  solve  the  x-momentum  and  energy  equations;  m^V'was  deter¬ 
mined  from  the  surface  mass  flux  which  was  specified,  and  ig  was 
determined  in  all  cases  in  a  manner  similar  to  that  used  in  Ref.  1  for 
laminar  flow,  namely  by  considering  the  x-momentum  equation  near  the 
outer  edge  of  the  boundary  layer.  The  value  so  obtained  was  adjusted 
to  prevent  large  variations  in  mg '  from  step  to  step  in  x,  and,  since 
mg'  controls  the  extent  of  the  region  over  which  the  boundary -layer 
equations  are  solved,  it  was  also  adjusted  to  assure  that  the  boundary- 
layer  equations  were  not  being  solved  over  a  region  of  essentially  in- 
viscid  flow. 

For  the  present  application,  the  primary  advantage  of  the  u  coordi¬ 
nate,  as  opposed  to  retaining  the  physical  coordinate  y,  is  that  the 
boundary  layer  always  lies  in  the  region  0  <  u  <  1,  while  in  terms  of  y 
the  extent  of  the  region  of  interest  is  unknown  and  generally  increases 
with  x.  If  the  variation  of  p  with  y  were  not  being  considered,  the  need 
for  a  knowledge  of  v  would  vanish,  and  neither  the  stream  function 
equation  nor  the  continuity  equation  would  need  to  be  solved. 


2.3  STAGNATION  POINT  EQUATIONS 

In  order  to  begin  the  numerical  solution  of  the  equations  represented 
by  Eq.  (17),  it  is  necessary  to  have  initial  profiles  of  u  and  H  as  func¬ 
tions  of  u.  These  can  be  provided  from  experimental  data,  approxima¬ 
tions,  or  other  sources,  as  in  Refs.  2  and  6  for  instance.  For  the  case 
of  flow  over  a  blunt-nosed  body,  however,  the  governing  equations  can 
be  reduced  to  a  two -point  boundary  -value  problem  at  the  stagnation 
point,  the  solution  of  which  can  provide  the  initial  profiles  needed  to 
start  the  general  solution. 

In  the  present  work,  the  stagnation  point  x-momentum  and  energy 
equations  were  formed  by  introducing  a  dependent  variable,  ?,  where 
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C  =  u/ue 


(22) 


and  then  considering  the  limiting  forms  which  the  equations  took  as  x 
approached  zero.  A  typical  term  in  this  treatment  is  the  following: 


lim 


x-*0  pei 


(23) 


where 


lim 

x-*0  r 

w 


In  these  equations,  (due/dx)0  is  the  stagnation  point  edge  velocity  grad¬ 
ient  which  was  evaluated  using  Euler's  equation  and  a  modified  Newtonian 
pressure  distribution. 


The  resulting  form  of  the  x-momentum  equation  at  the  stagnation 
point  is 


where 
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and 


T' 

W 


T; 


(i  + 


j  + 1 


o 


T'  =• 


(1 


y/rn)J 


(26) 


In  obtaining  du  the  y-momentum  equation  was  converted  to  x-u  form 
and  used  in  treating  those  terms  in  du  which  included  (9p/9u)x  as  a  coeffi¬ 
cient. 


Before  the  limiting  form  of  the  energy  equation  could  be  developed, 
it  was  necessary  to  multiply  the  entire  equation  by  ue.  The  resulting 
form  of  the  energy  equation  at  the  stagnation  point  is 

+  k  ['"(s)Jx  +  3ll  -  0  (27) 

where 

a,.  =  a 
11  u 


C[,  =  Cu/Pr 

3h  =  0  (28) 

Application  of  l'Hospital's  rule  to  the  definition  of  (j,  together  with 
the  introduction  of  the  definition  of  0  yields  the  following  equation  for  v 
at  the  stagnation  point 


v 


_ l _ 

( 1  +  Ky  )p{  1  +  y/rjJ 


+  y. 


'rn)J  nC  ~ 


K]  + 


(29) 
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Since  u  =  0  at  the  stagnation  point,  the  use  of  Eq.  (3)  yields  a  con¬ 
stant  static  pressure  at  that  location,  and  at  the  stagnation  point  y  can 
be  determined  from 


>  ■  /  izr 


(30) 


The  quantity  mg',  called  the  entrainment  rate  because  of  its  relation 
to  the  growth  of  the  stream  function  due  to  the  entrainment  of  fluid  into 
the  boundary  layer,  was  evaluated  at  the  stagnation  point  by  examining  a 
number  of  stagnation  point  boundary -layer  solutions  and  determining  that 


-30PcJIM\/Re 


yielded  a  value  which  was  generally  satisfactory. 


(31) 


2.4  FINITE-DIFFERENCE  EQUATIONS 

In  order  to  solve  the  coupled  partial  differential  equations  for  u  and 
H  represented  by  Eq.  (17),  this  equation  was  replaced  by  a  set  of  con¬ 
sistent  linearized  algebraic  equations.  The  set  is  of  tridiagonal  form 
and  was  solved  by  means  of  a  standard  algorithm  which  is  especially 
efficient  for  such  forms. 

As  mentioned  previously,  in  order  to  solve  the  finite -difference 
form  of  the  governing  equations  at  a  general  x  location,  it  was  necessary 
to  have  a  complete  solution  at  a  previous  x  station.  To  start  the  process 
this  initial  solution  was  obtained  from  experimental  data,  approximations, 
or  a  finite -difference  solution  of  the  stagnation  point  equations.  This  is 
indicated  schematically  in  Fig.  2  where  x  location  1  represents  a  loca¬ 
tion  where  the  solution  is  known  and  location  2  the  location  where  the 
solution  is  to  be  obtained. 

With  Eq.  (17)  in  the  form 

r  i  (dct\  _  af  _ i_ (ds\  0 

+  [_c,\df,>/x  bf<:fJ\<Wx  bf«f  +  bfcf  (32) 

the  following  difference  forms  obtained  from  Taylor  series  truncation 
were  employed. 
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*1,2  =  *1,1  +  °^x) 


Figure  2.  Typical  section  of  finite-difference  grid. 

Application  of  Eq.  (33)  to  Eq.  (32)  yields  the  following  general  form 
for  the  corresponding  finite -difference  equation 

*1,2  =  Af,I*I+l,2  +  Bf,  1*1-1, 2  *  cf,i  (34) 

where 

A  1  (  2  ,  Pf’*  \ 

Pf>I  +  A<a  j")  Aa  Awf  +  Acjj  J 

B  =  J_  /  2 _ PM  \ 

U  Df,l  yAeuf  +  Ao>pAw7  Aty|  +  ) 
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In  Eq.  (35)  it  is  indicated  that  af,  bf,  cf,  and  df  are  evaluated  at  x 
location  1,  and  this  proved  to  be  a  valid  and  effective  treatment  for 
laminar  flow  problems.  For  turbulent  flow  problems,  however,  diffi¬ 
culty  was  encountered  in  the  solution.  This  appeared  as  an  oscillation 
in  the  turbulent  viscosity  and  shearing  stress  from  one  x  location  to  the 
next,  in  the  neighborhood  of  the  changeover  from  the  inner-layer  vis¬ 
cosity  law  to  the  outer-layer  law.  To  alleviate  this  problem,  the  value 
of  the  viscosity  at  x  location  2  was  used  in  the  cf  terms.  Since  this  vis¬ 
cosity  is  not  known  beforehand,  this  treatment  necessitates  an  iteration 
for  the  solution  at  location  2. 

For  f  equal  to  each  of  u  and  H  an  equation  of  the  form  of  Eq.  (34) 
can  be  written  for  1  <  I  <  N,  where  N  is  the  number  of  grid  points  from 
the  surface  to  the  outer  edge,  including  both  end  points.  With  the 
boundary  conditions  as  described  previously,  the  equations  represented 
by  Eq.  (34)  were  solved  using  the  standard  algorithm  for  linear  tridiago¬ 
nal  systems,  as  described  in  Ref.  1.  After  the  values  of  u  and  H  were 
determined  at  x  location  2,  the  thermodynamic  and  transport  properties, 
p(y),  and  v(y)  were  determined.  Since  iteration  was  required,  for  the 
reason  previously  discussed,  this  was  done  until  convergence  was 
attained,  and  then  the  whole  procedure  was  repeated  at  the  next  x  loca¬ 
tion,  with  location  1  values  having  been  replaced  by  location  2  values. 

It  is  noteworthy  that  the  finite -difference  grid  used  required  neither 
Ax  nor  Au  to  be  a  constant;  in  particular,  an  w  grid  was  used  which 
caused  the  point  spacing  to  be  much  closer  near  the  surface  than  in  the 
outer  part  of  the  boundary  layer. 

At  the  stagnation  point  of  a  blunt  body  an  iterative  finite -difference 
method  was  used,  employing  difference  forms  for  the  derivatives  with 
respect  to  u  as  given  in  Eq.  (33).  All  data  except  the  values  of  the 
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dependent  variables  in  the  derivatives  with  respect  to  w  were  taken  from 
initial  approximation  or  previous  iteration  values  and  were  recomputed 
after  new  values  for  $  and  H  were  found.  The  iteration  was  then  re¬ 
peated  until  convergence  was  attained. 


3.0  RESULTS  OF  CALCULATIONS 


The  general  behavior  of  the  solution  method  described  in  the  pre¬ 
vious  section  was  investigated  by  performing  a  series  of  calculations  to 
investigate  grid  spacing,  convergence  criteria,  and  other  factors  in¬ 
volved  in  applying  the  finite -difference  technique.  It  was  determined 
that  a  practical  value  for  the  number  of  u  grid  lines  was  50,  and  the 
various  iterations  involved  were  deemed  to  be  converged  when  succes¬ 
sive  values  of  the  variables  involved  differed  no  more  than  0.  05  percent. 

All  computations  performed  in  this  investigation  were  done  in  double 
precision  on  an  IBM  Model  370/155  digital  computer,  using  approxi¬ 
mately  180,000  bytes  of  core  and  15  minutes  of  central  processor  time 
per  case. 


3.1  22.5-DEG  HYPERBOLOID  IN  MACH  10  FLOW 

One  test  of  the  present  method  was  to  make  comparisons  with  results 
obtained  in  the  investigation  reported  in  Ref.  4.  The  situation  was  that 
of  the  flow  of  air  over  a  0.  1-ft  nose-radius,  22.  5 -deg,  asymptotic  half¬ 
angle  hyperboloid.  The  free-stream  Mach  number  was  10,  and  the  free- 
stream  Reynolds  number  based  on  nose  radius  was  4,  000.  .  The  wall-to- 
free- stream  stagnation  temperature  ratio  was  0.2,  and  the  boundary - 
layer  flow  was  assumed  to  be  laminar.  For  this  case  the  boundary 
layer  occupies  about  one -third  of  the  region  between  the  body  and  the 
bow  shock,  and  the  effects  of  entropy  layer  swallowing  were  shown  in 
Ref.  4  to  be  significant. 

Figures  3,  4,  and  5  show  comparisons  of  results  from  the  present 
method  with  results  from  the  fully  viscous  shock  layer  method  of  Ref.  12. 
For  this  problem,  the  shock  layer  method  should  be  the  best  of  all  those 
which  are  available;  however,  shock  layer  methods  are  currently  quite 
restricted  in  their  range  of  applicability. 
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Figure  3.  Stanton  number  and  skin  friction  coefficient  distributions 
on  a  22.5-deg  hyperboloid  in  a  Mach  10  flow. 


- Present  Boundary-Layer  Results  with  Entropy  Swallowing,  dp/dy  *  0 

o  viscous  Shock  Layer  Method,  Ref.  12 


u/Um  H/H0 


Figure  4.  Velocity  and  total  enthalpy  profiles  at  x/rn  =  15 
on  a  22.5-deg  hyperboloid  in  a  Mach  10  flow. 
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- Present  Boundary-Layer  Results  with  Entropy  Swallowing,  dp/dy  t  0 

o  Viscous  Shock  Layer  Method,  Ref.  12 
□  Method  of  Characteristics,  Ref.  13 


Figure  5.  Static  pressure  and  normal  velocity  profiles  at  x/r„  =  15 
on  a  22.5-deg  hyperboloid  in  a  Mach  10  flow. 

Figure  3  compares  Stanton  number  and  skin  friction  coefficient  dis¬ 
tributions  from  the  present  method  and  the  viscous  shock  layer,  and 
excellent  agreement  between  the  two  methods  is  exhibited. 

Figure  4  shows  velocity  and  total  enthalpy  profile  data  from  the  two 
methods  at  a  surface  location  of  x/rn  =  15.  The  boundary -layer  solution 
was  terminated  at  a  y  location  such  that  B/B^  had  attained  a  value  close 
to  unity  at  y/ye  =  0.  9;  the  shock  layer  data  extend  out  to  the  computed 
location  of  the  bow  shock  wave.  These  profile  data  also  show  good 
agreement  between  the  two  methods. 

Figure  5  compares  static  pressure  and  normal  velocity  data  from 
the  present  method  and  the  shock  layer  method  at  x/rn  =15.  In  addition, 
pressure  data  obtained  using  the  inviscid  method  of  characteristics  treat¬ 
ment  of  Ref.  13  is  shown.  There  is  reasonable  agreement,  considering 
the  expanded  abscissa  scale,  among  the  three  methods  where  they  over¬ 
lap;  however,  noting  that  Eq.  (3)  requires  9p/9y  >  0  for  this  case,  it  is 
obvious  that  the  outer  region  of  decreasing  static  pressure  could  not  be 
obtained  using  this  model  for  treating  the  flow  in  such  a  region. 

Obviously,  as  mentioned  in  Ref.  4,  the  static  pressure  variation  shown 
here  is  negligibly  small;  it  has  been  considered  only  for  demonstrating 
the  capability  for  treating  it.  The  normal  velocity  data  show  good  agree¬ 
ment  between  the  present  results  and  those  obtained  by  the  method  of 
Ref.  12. 
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The  results  discussed  in  this  subsection  serve  to  validate  the  basic 
method  which  has  been  presented  for  treating  laminar  boundary -layer 
flows,  and  to  demonstrate  the  capability  (and  some  limitations)  for  com 
puting  the  static  pressure  variation  with  y  and  the  normal  component  of 
velocity  in  boundary -layer  flows. 


3.2  CONE-OGIVE-CYLINDER  IN  MACH  7.2  FLOW 

In  Refs.  14  and  15,  Horstman  and  Owen  present  detailed  measure¬ 
ments  of  the  flow  in  the  turbulent  boundary  layer  on  a  slender  cone¬ 
ogive-cylinder  body  in  a  Mach  7.  2  flow.  They  report  measurements  of 
surface  pressure,  skin  friction,  and  heating  rate,  profiles  of  various 
mean  flow  parameters,  and  results  of  hot-wire  probe  measurements. 
Figures  6  and  7  and  Table  1  compare  results  from  the  present  investi¬ 
gation  with  the  data  obtained  by  Hortsman  and  Owen  for  all  except  their 
hot-wire  results. 


b.  Body  and  bow  shock 

Figure  6.  Surface  pressure  distribution,  body,  and  bow  shock  for 
'  cone-ogive-cylinder  in  a  Mach  7.2  flow. 
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Table  1.  Boundary- Layer  Parameters  for  Cone-Ogive-Cylinder  in  a 
Mach  7.2  Flow 


Measured 

Computed 

Measured 

Computed 

Measured 

Computed 

Station 

1 

1 

2 

2 

3 

3 

x,  ft 

3.80 

3.80 

5.  80 

5.80 

7.80 

7.80 

4,  ft 

0.0558 

0.053 

0.  082 

0.079 

0. 1083 

0.  104 

4*,  ft 

0.  0284 

0.0287 

0.0427 

0.0432 

0.0531 

0.  0557 

8,  ft 

0.  00176 

0.00174 

0.00273 

0.00271 

0.00390 

0.00360 

cfe  x  103 

0.90 

0.884 

0.  855 

0.  854 

0.801 

0.835 

5te  x  103 

0.510 

0.612 

0.458 

0.569 

0.484 

0.  543 

Notes:  1.  4,0  include  transverse  curvature 


2.  4  -  y  at  Po/(Po>y^„  *  00 

3.  cfc  =  2tw/p6u^ 

4.  Ste  !  <1  w /  [Peuc^*  9  Hg  -  H^) ] 

The  body  was  approximately  12  ft  long,  being  basically  a  4-in.-diam 
cylinder  capped  by  a  10-deg  half-angle  cone  followed  by  short  ogive  sec¬ 
tion;  the  cylinder  began  approximately  2  ft  past  the  cone  apex.  The  flow 
underwent  transition  from  laminar  to  turbulent  flow  beginning  at  x  =  1.2 
ft  and  ending  at  x  =  2.  6  ft,  and  the  surface-to-free-stream  total- 
temperature  ratio  was  0,465. 

Figure  6  shows  the  surface  pressure  distribution,  the  body  shape, 
and  the  bow  shock  shape  for  the  case  of  Horstman  and  Owen.  The  sur¬ 
face  pressure  distribution  computed  using  the  method  of  characteristics 
technique  of  Ref.  13  agrees  reasonably  well  with  the  experimental  data, 
although  the  computed  value  is  about  10  percent  higher  than  that  mea¬ 
sured  at  the  last  measuring  point.  In  Fig.  6b,  the  computed  boundary- 
layer  thickness  is  seen  to  be  less  than  10  percent  of  the  distance  from 
the  body  to  the  bow  shock  at  the  locations  where  bound  ary -layer  profile 
and  surface  heating  rate  and  skin  friction  measurements  were  made. 

Also  shown  in  Fig.  6b  is  the  location  (at  each  survey  station)  out  to 
which  the  computed  static  pressure  was  essentially  constant.  Because 
the  body  longitudinal  curvature  is  zero  over  the  cylinder,  the  boundary- 
layer  computations  assumed  9p/9y  =  0  in  that  region.  Inclusion  of  longi¬ 
tudinal  curvature  effects  over  the  ogive  section  had  no  significant  effect 
on  the  flow  at  the  survey  stations. 

Table  1  compares  boundary -layer  thickness  data,  integral  quanti¬ 
ties,  and  surface  heating  and  friction  parameters  as  measured  with 
those  computed.  The  computed  and  experimentally  determined  thickness 
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and  integral  parameters  are  in  excellent  agreement.  The  computed 
coefficient  of  skin  friction  is  within  ±4  percent  of  the  experimental  data 
at  each  station.  The  computed  Stanton  number  is  approximately  20  per¬ 
cent  greater  than  the  experimental  data;  however,  the  experimental 
data  do  not  show  the  expected  monotonically  decreasing  behavior. 


Figure  7  shows  boundary-layer  profiles  at  survey  station  2.  The 
comparisons  between  the  computed  results  and  the  experimental  data 
were  very  similar  to  these  at  all  three  survey  stations,  and  these  are 
presented  as  typical.  There  is  good  agreement  between  computation 
and  experiment  in  the  case  of  the  velocity,  total  enthalpy,  and  Mach 
number  profiles.  The  density  profiles  agree  well  except  near  the  outer 
edge  of  the  boundary  layer,  where  the  computed  data  reach  a  value 
almost  15  percent  greater  than  that  measured. 


3.3  MACH  8  NOZZLE  FLOW 

In  this  section  are  presented  comparisons  between  experimental 
data  and  computations  for  the  case  of  supersonic  flow  in  a  nozzle.  The 
experimental ^data  were  taken  in  the  Arnold  Engineering  Development 
Center  vonKarman  Facility's  Hypersonic  Wind  Tunnel  B  at  an  exit  Mach 
number  of  approximately  8.  The  axisymmetric  contoured  nozzle  is 
approximately  22  ft  long,  with  a  throat  radius  of  1.  60  in.  and  an  exit 
radius  of  about  2  ft.  Tabulated  geometry  data  are  given  in  Ref.  16. 

The  test  condition  considered  was  one  for  which  the  total  pressure 
was  605  psia  and  the  total  temperature  132 1°R,  At  the  survey  station 
where  the  experimental  data  were  taken,  the  unit  Reynolds  number  out¬ 
side  the  boundary  layer  was  2.  69  x  10^/ft,  and  the  corresponding  Mach 
number  was  8.  02.  These  data  and  the  experimental  data  which  will  be 
presented  were  taken  from  Ref.  17,  Surface  temperature  data  used  in 
the  computations  were  taken  from  estimates  given  in  Ref.  18. 

The  inviscid  data  which  provided  the  boundary  conditions  used  along 
the  outer  edge  of  the  boundary  layer  were  determined  by  an  isentropic 
expansion  from  reservoir  conditions  to  a  known  pressure  at  a  given 
value  of  x.  The  pressure  distribution  (in  x)  was  obtained  from  one¬ 
dimensional  stream  tube  relationships.  Viscous  interaction  effects 
were  considered  by  forming  an  effective  area  distribution  which  took 
into  account  the  bound  ary -layer  displacement  thickness.  This  required 
an  iteration  with  respect  to  the  displacement  thickness,  in  a  manner 
similar  to  that  described  in  Ref.  19. 
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Figure  8  shows  the  results  of  the  iteration  procedure,  beginning 
with  an  initial  approximation  to  the  displacement  thickness  distribution 
which  was  a  bit  high,  and  converging  reasonably  well  in  2  iterations. 
Because  of  the  difficulties  in  defining  and  computing  a  displacement 
(or  momentum)  thickness  when  considering  longitudinal  curvature  and 
normal  pressure  gradient  effects,  these  effects  were  neglected  in 
these  nozzle  bound  ary -layer  computations.  (The  initial  boundary- 
layer  solution  was  performed  with  and  without  the  longitudinal  curva¬ 
ture/normal  pressure  gradient  effects,  and  the  net  effect  on  the  solu¬ 
tion  was  negligible.) 


Figure  8.  Displacement  thickness  iterations  for  Mach  8  nozzle  case. 

Table  2  lists  the  measured  and  computed  displacement  and  momen¬ 
tum  thickness  values  at  the  survey  station.  The  computed  displacement 
thickness  is  approximately  fifteen  percent  greater  than  the  experimen¬ 
tally  inferred  value,  and  the  computed  momentum  thickness  almost  thirty 
percent  less  than  the  experimental  value.  It  might  be  noted  that  the  re¬ 
sults  of  nozzle  bound  ary -layer  computations  similar  to  the  present  ones 
were  given  in  Ref.  20,  and  the  general  nature  of  the  agreement  between 
the  computed  and  measured  data  was  essentially  as  obtained  in  this  inves¬ 
tigation. 
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Table  2.  Boundary-Layer  Parameters  at  Exit  of  Mach  8  Nozzle 

g*,  ft  e,  ft 

Measured  0.  198  0.0221 

Computed  0.227  0.  0160 

Figure  9  shows  comparisons  between  computed  and  measured  im¬ 
pact  pressure  and  total  enthalpy  profiles  at  the  survey  station.  Reason¬ 
able  agreement  is  shown  between  the  computations  and  the  experimental 
data.  The  computed  profiles  extend  as  far  in  y  as  the  calculations  were 
carried,  while  the  experimental  data  actually  continued  beyond  the 
region  shown  but  with  constant  values  of  the  pitot  pressure  and  total 
enthalpy. 


Figure  9.  Impact  pressure  and  total  enthalpy  profiles  at 
exit  of  Mach  8  nozzle. 


3.4  RAMP  IN  MACH  3.54  FLOW 

In  Ref.  21,  Sturek  presents  the  results  of  an  experimental  investi¬ 
gation  of  the  turbulent  boundary-layer  flow  over  a  two-dimensional  isen- 
tropic  compression  ramp  with  significant  longitudinal  curvature  and 
normal  pressure  gradient  effects.  The  ramp  was  affixed  to  the  floor  of 
the  test  section  of  a  supersonic  wind  tunnel  adjusted  to  give  a  Mach  3.  54 
uniform  flow  approaching  the  ramp.  The  condition  considered  in  this 
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study  was  one  for  which  the  total  temperature  was  560°R  and  the  total 
pressure  was  48.  34  psia.  The  ramp  surface  was  essentially  adiabatic, 
and  the  test  section  boundary  layer  approaching  the  ramp  was  approxi¬ 
mately  1  in.  thick. 

In  order  to  begin  the  computations  one-seventh  power  profiles  of 
u/ue  and  (H-Hw)/(He-Hw)  were  assumed,  giving  the  comparison  shown 
in  Fig.  10  between  the  assumed  and  measured  velocity  and  total  enthalpy 
profiles  approaching  the  ramp.  The  velocity  profiles  agree  quite  well, 
and  the  total  enthalpy  profiles  agree  to  within  about  2  percent,  although 
the  assumed  profile  does  not  exhibit  the  overshoot  to  be  expected  in  the 
outer  region  of  the  boundary  layer  on  an  adiabatic  surface.  At  the 
location  where  the  computations  were  begun,  the  static  pressure  was 
constant  across  the  boundary  layer. 


Figure  10.  Assumed  and  measured  profiles  near  the  beginning  of  a  ramp  in 
a  Mach  3.54  flow. 


The  measured  pressure  distribution  along  the  ramp  was  used  to  pro¬ 
vide  the  necessary  data  for  determining  the  inviscid  flow  over  the  ramp, 
according  to  the  scheme  discussed  in  Section  2.  1.  The  coordinates  of  the 
ramp  were  fitted  by  a  curve  of  the  form  y  =  ax*3  (where  x  and  y  are  rec¬ 
tangular  cartesian  coordinates),  and  the  corresponding  longitudinal  radius 
of  curvature  was  determined  analytically. 
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Figure  11  shows  computed  and  measured  profiles  of  velocity,  total 
enthalpy,  and  static  pressure  at  a  survey  station  approximately  8  in. 
back  on  the  ramp. 


-Computed 


y. 


Figure  11.  Computed  and  measured  velocity,  total  enthalpy,  and 
static  pressure  profiles  at  survey  station  on  ramp 
in  a  Mach  3.54  flow. 


Good  agreement  is  shown  between  the  measured  and  computed  data, 
except  that  the  computed  total  enthalpy  profile  still  shows  no  overshoot. 
In  particular,  the  static  pressure  variation  is  felt  to  be  quite  well  repre 
sented  by  the  model  used  in  the  computations. 


Figure  12  compares  computed  and  measured  surface  shearing  stress 
distributions  along  the  ramp.  The  two  sets  of  data  agree  reasonably  well, 
but  the  experimental  data  increase  less  rapidly  over  the  forward  portion 
of  the  ramp  and  more  rapidly  over  the  aft  portion  of  the  ramp  than  the 
computed  shearing  stress. 
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Figure  12.  Computed  and  measured  surface  shearing  stress  distribution 
on  ramp  in  a  Mach  3.54  flow. 


4.0  SUMMARY  AND  CONCLUSIONS 


The  compressible  laminar  and  turbulent  bound  ary -layer  equations 
for  two-dimensional  and  axisymmetric  flow  have  been  presented  in  a 
form  which  includes  transverse  curvature,  longitudinal  curvature,  and 
normal  pressure  gradient  effects.  The  equations  were  first  presented 
in  physical  surface  coordinates  and  then  transformed  using  the  u2 
transformation  to  a  form  for  which  a  finite -difference  solution  procedure 
was  developed.  The  stagnation  point  form  of  the  boundary-layer  equa¬ 
tions  in  x-u  coordinates  was  also  presented,  to  the  author's  knowledge 
the  first  time  this  has  been  done. 

The  method  presented  for  treating  boundary -layer  flows  was  tested 
by  comparing  results  from  the  present  method  with  both  experimental 
data  and  results  from  other  methods  of  calculation.  Generally  good 
agreement  was  obtained  for  a  variety  of  cases,  including  both  internal 
and  external  boundary  layers  and  a  case  in  which  normal  pressure 
gradient  and  longitudinal  curvature  effects  were  significant. 

The  method  presented  herein  would  be  difficult  to  apply  to  situations 
in  which  boundary -layer  displacement  influences  needed  to  be  considered. 
This  is  because  of  the  difficulty  which  would  be  encountered  in  defining 
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and  computing  a  displacement  thickness  when  normal  pressure  gradient 
and  longitudinal  curvature  effects  were  included.  The  best  approach  for 
such  cases  would  be  the  development  of  a  fully  viscous  shock  layer 
method  of  more  general  applicability  than  that  given  in  Ref.  12. 
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NOMENCLATURE 

Terms  in  finite -difference  equations 
Mixing  length  constant 

Terms  in  general  form  of  governing  equations 

.  2 

Edge  skin  friction  coefficient,  2Tw/peue 

2 

Free-stream  skin  friction  coefficient,  2TW/PooUa, 

General  variable,  especially  u  or  H 

Stagnation  enthalpy 

Mixing  length  constant 

Mach  number 

Mass  flux  term 

Number  of  cross-stream  grid  points 

tfe  “  '^w 
Pressure  1 

Pitot  pressure 

Prandtl  number 

Heat -transfer  rate 

Longitudinal  radius  of  curvature 

Free-stream  Reynolds  number  based  on  nose  radius 

Radial  coordinate 

Nose  radius  of  curvature 

Edge  Stanton  number 

Free-stream  Stanton  number,  q-w/Pootla^Ha,  -  Hw) 
Temperature 

Terms  in  stagnation  point  equations 
Free-stream  velocity 
Velocity  in  x  direction 
Velocity  in  y  direction 
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X 

y 

Ym 

z 
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6 

6* 
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#c 
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Coordinate  along  the  surface 
Coordinate  normal  to  the  surface 
Value  of  y  at  u/ue  =  0.  99 
Axial  distance 
Intermittency  factor 
Boundary -layer  thickness 
Boundary-layer  displacement  thickness 
u/ue 

Momentum  thickness 
Longitudinal  curvature,  1/R 
Mixing  length  constant 
Viscosity 
Mass  density 
Shearing  stress 
Surface  slope 
Stream  function 
Transformed  stream  function 


SUBSCRIPTS 

e  Outer  edge  of  boundary  layer 

i  Laminar 

t  Turbulent 

w  Wall  or  surface 

o  Stagnation  conditions 

-  Free-stream,  or  reservoir  for  internal  flow 


34 


